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We consider relaxation of a single electron spin in a nanotube quantum dot due to its coupling 
to flexural phonon modes, and identify a new spin-orbit mediated coupling between the nanotube 
deflection and the electron spin. This mechanism dominates other spin relaxation mechanisms in 
the limit of small energy transfers. Due to the quadratic dispersion law of long wavelength flexons, 
a; oc g^, the density of states dq/doj oc uj"^^'^ diverges as a; — ^ 0. Furthermore, because here the spin 
couples directly to the nanotube deflection, there is an additional enhancement by a factor of l/g 
compared to the deformation potential coupling mechanism. We show that the deflection coupling 
robustly gives rise to a minimum in the magnetic field dependence of the spin lifetime Ti near an 
avoided crossing between spin-orbit split levels in both the high and low-temperature limits. This 
provides a mechanism that supports the identification of the observed Ti minimum with an avoided 
crossing in the single particle spectrum by Churchill et al.[Phys. Rev. Lett. 102, 166802 (2009)]. 

PACS numbers: 72.25.Rb, 71.70.Ej, 85.35.Kt 



I. INTRODUCTION 

Due to their outstanding mechanical properties and 
versatile electrical characteristics, carbon nanotubes of- 
fer an exciting platform both for studies of fundamental 
physical phenomena and for a variety of potential appli- 
cations. The relatively small nuclear charge of carbon 
and the low natural abundance of carbon isotopes with 
nonzero nuclear spin suggest that the spin-orbit and hy- 
perfine interactions, which are the main sources of elec- 
tron spin relaxation in GaAsi— , should be weak in car- 
bon nanotubes. Thus in recent years the electronic spin 
properties of nanotubes have gained wide attention for 
potential applications in spintronics and quantum com- 
puting. Furthermore, the availability of isotopically pu- 
rified starting materials opens the possibility of growing 
^^C (nuclear spin / = 0) and ^■^C (nuclear spin / ~ 1/2) 
nanotubes to study the behavior of electron spins in the 
presence or absence of a nuclear spin bathSi. 

As techniques for preparing ultra-clean samples and 
studying them in cryogenic environments have become 
available, few-electron quantum dots have emerged as a 
powerful tool for studying the electron spin properties 
of nanotubes^"—. In the experiment by Kuemmeth et 
dXA, single electron and single hole quantum dot spec- 
tra were shown to display the characteristics of coher- 
ent coupling between the electron's spin and its orbital 
magnetic moment^^— , see FigUJi. Experiments in the 
spin-blockade regime of few electron double quantum 
dots^ d'^d^ have also yielded important information about 
spin relaxation in nanotubes. In particular, the experi- 
ment by Churchill et ali^ demonstrated a minimum of 
the spin lifetime Ti near the narrow avoided crossing be- 
tween levels circled in Fig|TJ Below we will identify a 
mechanism of spin relaxation in nanotube quantum dots 
that is based on the coupling of an electron's spin to the 
nanotube's deflection. This mechanism provides a deeper 
understanding of the Ti minimum. Our theory is devel- 
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FIG. 1: Spin relaxation due to defiection-coupling to bend- 
ing mode phonons. a) Single electron quantum dot energy 
level diagram with the parameter values of Ref . 0] . The four- 
dimensional subspace is spanned by the states | t s ) labeled 
by the valley index r = ±1 and the spin projection s = ±1 
along the laboratory z-axis, see Eqs.© and Q. For the mag- 
netic field B = 0, the upper (lower) Kramers doublet includes 
states I 1, 1 > and | -1, -1) (| 1, -1 ) and | -1, 1 )). Due to the 
flexon density of states singularity as a; — > 0, we focus on spin 
relaxation rates between levels with small energy splittings 
in the lower Kramers doublet (denoted by E±), and in the 
narrow upper avoided crossing (dashed circle), b) Gedanken 
experiment to illustrate the coupling mechanism. In the ab- 
sence of spin-orbit coupling, the electron spin remains fixed 
in the laboratory frame irrespective of the nanotube's mo- 
tion. When the strength of spin-orbit coupling Aso is much 
greater than the rate of motion fl, coupling between the spin 
and orbital moments keeps the electron spin aligned with the 
direction of the tube axis, c) For a nanotube with constrained 
ends, bending mode phonons cause spatial variations in the 
direction of the nanotube axis t{z), which locally couple to 
the electron spin as described above. 

oped for a single electron quantum dot, with the corre- 
sponding energy level diagram shown in FiglTJi. While 
we envision that the basic results are more general, a full 
investigation of the details for many-electron systems is 
beyond the scope of this paper. 

In order for spin relaxation to occur, energy must be 
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transferred from the spin to the environment, which for 
confined electrons typicahy involves the phonon bath of 
the host crystal. Due to strong coupling, phonons play 
an especially important role in nanotubes^^^— . Crucially 
for us, clean nanotubes possess quadratically dispersing 
bending mode o^^'-'^^i^^'^" with uj{q) oc q^. The corre- 
sponding density of states for each such mode has a Van 
Hove singularity dq/du cx '^j \fi^ at low energies (w — > 0). 
Thus any relaxation mechanism that relies on coupling 
to such phonons is expected to become especially efficient 
for small energy transfers. 

Indeed, Bulaev et alJ^ considered spin-orbit coupling 
to the bending mode deformation potential, and used the 
singularity in the density of states to predict a spectac- 
ular enhancement of 1/Ti oc l/\/uj near the upper anti- 
crossing point of Fig[T^ in the high temperature regime. 
In the low temperature regime, however, this mechanism 
predicts a suppression of spin relaxation 1 /Ti cx ^/uj. Be- 
cause the experiments of Ref.jH were performed in an 
intermediate regime huj/ks ~ T ^ 100 mK, it is dif- 
ficult to confidently link the Ti minimum to the flexon 
deformation potential coupling. The deflection-coupling 
mechanism described below resolves this ambiguity, as it 
ensures a Ti minimum for all temperatures. 

Borysenko et ali^i also invoked the density of states 
singularity, but rather than the deformation potential 
they considered spin-flcxon coupling through the Zeeman 
interaction due to the g-factor anisotropy Aq. This cou- 
pling is weak because /S.g sa 10~^, see Ref.[22[ and foot- 
note [131 • Moreover, because the orbital moment^ /iorb 
was not taken into account, the upper avoided crossing 
did not appear in their spectrum and the theory was only 
applied to the low-magnetic field Kramers doublets. 

Here we consider a different coupling mechanism be- 
tween flexural modes and electron spin. Similar to the 
mechanism of Borysenko et ali^, the coupling is pro- 
portional to the nanotube deflection, arising from local 
changes of the direction of the nanotube axis in the global 
(laboratory) reference frame. Unlike their mechanism, 
however, the spin-orbit interaction of Refs.flQ is intrin- 
sically sufflcient to provide the coupling, which thus op- 
erates even in zero external magnetic fleld. Therefore, 
our mechanism is efficient near both the upper avoided 
crossing and the narrow Kramers doublets, is not limited 
by small Ag, and is parametrically stronger than the de- 
formation potential mechanism at low energies (small q) 
because it is proportional to the deflection 6{z) rather 
than the deformation e ~ 86 /dz ^ qO. 

The physical origin of the coupling mechanism can 
most easily be understood through a thought experiment, 
illustrated in FigUJ). Imagine rigidly twirling a nanotube 
in free space about an axis perpendicular to the tube axis. 
In the absence of spin-orbit coupling, an electron spin ini- 
tially oriented in some direction in the global (laboratory) 
reference frame will remain polarized in the same direc- 
tion as the nanotube rotates around it. Now consider the 
opposite case of strong coupling between the electron spin 
and orbital magnetic moments. If the electron is initial- 



ized to a state with the orbital and spin moments aligned 
along the tube axis, then as the tube slowly rotates, the 
electron spin will adiabatically follow the direction of the 
tube axis. For weaker coupling and/or in the presence 
of an external magnetic field which provides a preferred 
direction for the electron spin in the laboratory frame, 
non-trivial spin-dynamics will be generated. 

In the more relevant situation of a freely vibrating, 
electrically-contacted, suspended or substrate-supported 
nanotube with immobilized ends, bending mode phonons 
can accomplish the same effect locally, section-by-section 
throughout the tube, as indicated in FiglTJ;. Below we 
treat this effect perturbatively for small deflections. To 
make connection to the Gedanken experiment depicted 
in FigHJ; where the spin follows the direction of the tube 
axis, we note that the clcctron-flexon coupling Hamil- 
tonian -ffg-ph [sec Eq.(l5l) and accompanying discussion 
below] used throughout this paper depends only on the 
instantaneous deflection of the nanotube. Because the 
local coupling strength depends on the deflection angle, 
it is proportional to the first derivative of the nanotube 
displacement along the tube axis z, and hence to the first 
power of the phonon momentum q. 

The plan for the rest of the paper is as follows. In 
section [III we present the model in detail and define the 
transition rates to be calculated. Then in section IIIII 
we calculate the spin relaxation rate between states of 
the field-split lower Kramers doublet for weak external 
field B K, Q. In section ITVl we calculate the spin relax- 
ation rate between the two nearly degenerate states at 
the upper avoided crossing of FiglT^ and demonstrate 
the scaling of the Ti minimum in the high-temperature 
and low-temperature limits. The main results of sections 
IIIII and HVl are summarized in Eqs.([T71) and respec- 
tively. Last, in section IVl we discuss how this picture is 
modified in "dirty" tubes where a substrate or coating 
may lead to flexon localization. 



II. SPIN PHONON COUPLING 

We consider a single electron conflncd in a semicon- 
ducting (narrow-gap or large-gap) nanotube quantum 
dot, which in the leading approximation is described by 
the Hamiltonian 



Hd = hvj- 



T^kcai - icr2 — 
dz 



V{z), 



(1) 



where vf ~ lO^m/s is the Fermi velocity of graphene, 
T3 is the isospin Pauli matrix with eigenvalue r in valley 
Kt, where Ki ~ 27r/3a(l, and K^i = — A'l arc the 
K and K' points of the graphene Brillouin zone and a is 
the lattice constant, cti and a2 are pseudospin (sublattice 
space) Pauli matrices, and V{z) is a conflning potential 
in the longitudinal direction z. The transverse momen- 
tum hkc, which sets the band gap Eg = 2hvpkc for free 
electrons when V{z) = 0, is determined by the boundary 
condition and warping effects that arise from rolling the 
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graphene sheet into a cyhnder. For large-gap nanotubes, 
fee is proportional to 1/i?, where R is the radius of the 
tube. For narrow-gap tubes, kc oc cos Bry/i?^ is non-zero 
due to curvature of the graphene shee1j2^. Here rj is the 
winding angle of the tube^ i^^i^^ . To keep the discussion 
general, for now we avoid imposing any particular form 
of the longitudinal confining potential. 

Due to the 2-fold real spin and 2-fold isospin symme- 
tries, all eigenstates of Hd are four-fold degenerate. In 
particular this applies to the ground state: 



Hd\T s)^Eo\t s] 



(2) 



where \t s) — \t)® \ s) is factorized into orbital and spin 
parts satisfying Tajr) = t|t), Sz\s) = s\s). Here is 
the real spin Pauli matrix associated with the z direction 
in the laboratory frame. Note that | t ) describes both 
the quantized longitudinal and transverse orbital motion 
of an electron in valley Kr- The specific expression for 
Eq depends on the structure of the nanotube and on the 
potential V{z). Below we set the zero of energy to be 

We now add another piece Hs to the full system Hamil- 
tonian H = Hd + Hs, that is central to this paper. It 
includes the coupling of an external magnetic field B to 
the electronic orbital and spin magnetic moments fiorb 
and ^B, spin-orbit coupling, and intervalley scattering, 

Hs = ^Y^T^ii ■ s) + Akk'Ti - (J-orbTsii ■ B) + fiB{s ■ B), 

(3) 

where t — (ii^,iy,iz) is the local tangent unit vector at 
each point along the tube, t — t{z), with components de- 
fined in the laboratory frame, and Aso a-nd Akk' are the 
spin-orbit and intervalley coupling matrix elements. The 
constant Aso in our model Hamiltonian Hs absorbs both 
currently known spin-orbit coupling terms^i^rii. Orbital 
angular moments were discovered by Minot et al. in 
Ref . Q , and spin-orbit coupling was found and measured 
by Kuemmeth et al. and Churchill et al. in Refs.jlQ. 
We assume that at rest the tube is straight and oriented 
along the laboratory z axis. 

Coupling to flexural phonons appears when we take 
into account the fact that the tangent vector is actually 
an operator depending on the phonon displacement co- 
ordinates u{z) at each point along the tube. For small 
amplitudes, long-wavelength deflections are described by 
t{z) = z -\- du±{z)/dz, where u±^ is the nanotube dis- 
placement perpendicular to z. Substituting this expres- 
sion for t into Eq.Q, we split the Hamiltonian Hs into 
a zero-order unperturbed part^ 



A: 



Ho = —^T^s^ + Akk'Ti - HorbTsB^ + ^b(s • JB), (4) 
and a spin-phonon coupling part 



H 



A, 



s — ph 



SO 



dz 



( du± , 

T3 - Morh I -J^ ■ B ]T3. (5) 



As usual in the theory of condensed matter systems, and 
in the spirit of the Born-Oppenheimer approximation, 
the spin-phonon coupling Hamiltonian Hg-ph depends 
only on the instantaneous deflection u±{z). In terms 
of the flexon creation and annihilation operators aj^, and 



du± 

dz 



(^) 



2pLzUJq 



(6) 



where a distinguishes two bending mode polarizations 
Xq, Lz is the nanotube length, q and ujq are flexon mo- 
mentum and frequency, and p is mass per unit length. As 
shown in Refs. [l5lll6ll2G |. uj{q) = Pq^ for small g; because 
of the nanotube symmetry, the two polarization modes 
are degenerate. 

The phonon normal mode profiles and frequencies de- 
pend in a non-universal way on the details of the nan- 
otube's environment (e.g. on boundary conditions, and 
on the presence of disorder or coating particles on the 
tube's surface). In Eq.® we capture the generic be- 
havior of the mechanism by assuming plane- wave normal 
modes. However, below we write the main results (|17p 
and (PD)) in terms of a form-factor M{q), Eq. ([5]), which 
can be generalized to other normal mode profiles. Indeed, 
a more general discussion of the flexon normal modes is 
provided in SecfVl 

Similar comments can be made regarding the influence 
of electron-electron interactions in many-electron quan- 
tum dots. A homogeneous background of closed electron 
shells should merely renormalize the coefficients with- 
out influencing the basic qualitative results. However, 
a strongly inliomogeneous electron background includ- 
ing charge puddles may influence electron-flexon coupling 
more profoundly. Such a regime is outside the scope of 
our investigation. 

The spin Hamiltonian Hs is defined on the four- 
dimensional subspace spanned by the eigenstates | t s ) 
of Hd, Eq.®. The orbital term porbTz{t ■ B) results from 
a shift of kc due to the vector potential associated with 
B. This renormalizes the gap Eg, which in principle 
affects the longitudinal confinement (see Appendix . 
Because this shift is small in the regime of interest for 
experiments p,orbii ■ B)/Eg < 1/30, we ignore its 

effect on the longitudinal motion. By ignoring coupling 
to higher orbital levels we omit terms which are small 
in the inverse level spacing, but in return obtain gen- 
eral analytical formulas which allow us to clearly extract 
the essential physics underlying the deflection coupling 
mechanism. 

Armed with the perturbation ([5]), we calculate transi- 
tion rates between eigenstates of the zero-order Hamilto- 
nian Hd + Hq using Fermi's Golden Rule: 



dq 



(7) 



nw=Ef-Ei 



Here | V'i ) and | ■0/ ) are the initial and final states satisfy- 
ing [Hd + Hq] \ipn) = En {ipn), and the over-bar indicates 
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averaging over the thermal phonon distribution. 

As discussed further in Appendix [XI due to the fact 
that the orbital states | r ) with r = ±1 are time-reversal 
conjugate, the density n(zo) = (t | S{z — zq) | t), which 
is a real scalar, is independent of the isospin index r. 
Using this fact and the property that Eq.® is diagonal in 
the isospin index r, any matrix element {ipf \ i?s-ph | V'l ) 
includes the longitudinal form factor 



M{q) = / dzn{z) e 



iqz 



(8) 



which depends on the specific form of the longitudinal 
confinement, but otherwise does not depend on the com- 
position of I "04 ) s-nd I "0/ ) in terms of the basis states 
{\rs)). 

In Appendix El we calculate M{q) for square- well con- 
finement. For small momentum g — >■ 0, the form factor 
M{q) — > 1, simply reflecting the wave function normal- 
ization. For large phonon momentum qL^ ^ 27r, where 
Ld is the electronic length of the quantum dot, the fac- 
tor e"?^ oscillates rapidly relative to the electron wave 
function and leads to a strong suppression of the matrix 
element, as will be discussed in section [V] Note that we 
distinguish between L,, the unconstrained length of the 



nanotube relevant for phonon normal modes, and Ld, 
the quantization length of the electronic wavefunction 
controlled by electrostatic gates, which in general are in- 
dependent quantities (see e.g. Ref.[3). 



III. RELAXATION RATE IN THE LOWER 
KRAMERS DOUBLET 



Due to the flexon density of states singularity as w — >■ 0, 
we focus our attention on spin relaxation between states 
with small energy splittings. We begin by considering 
relaxation between the lower Kramers pair of states (see 
FigH^) in a small applied magnetic field. Although time 
reversal symmetry prevents relaxation from occurring at 
B EE I J5| = (Van Vleck cancellation, see Refs.jHigSl]), 
the large density of states contributes to a steep rise of 
the relaxation rate as B is increased away from zero. 

To apply Eq.© for the relaxation rate, the first step is 
to find the two lowest energy eigenstates of Hamiltonian 
Hf). We take Aso > 0, as observed in Ref.Q- Diago- 
nalization is accomplished in two steps: we first apply a 
unitary transformation 



J 



U 



1 /A + A 



so 



2A 



1 - 



A.^ 



KK' 



SzT3 1 



\Akk' 



A- Aso 



2A 



I A,. 



KK' 



- iSzT2 1 



^KK' 



^KK' 



(9) 



to switch to a frame where UHq{B = 0)?/^ = —^^-'''3 ^® 

For Aso < 0, a 



Here A = V^Iq 



diagonal 

different unitary transformation must be chosen. 



4A^K' 



After the transformation, the zero-order Hamiltonian 
UHaW for B 7^ reads 



= - 



A Aso 

-yTaSz - ^J-orbBz—^Ts + HbB^Sz 



r. r, ^KK' jj AaAT'AsO „ n |AaA'| „>, 
-2fJ.BB± T2Sy- fiBBj_— —Sx-2^1orbBz T TiSz, (10) 



'^KK 



'lA 



r 



where without loss of generality we choose the perpendic- 
ular component of the magnetic field B±^ to be oriented 
along the x direction, B^ = B±,By = 0. 

To linear order in B/A, the lowest two energy eigen- 
values are E± = — (-j i with 



A 



(/Xor&AsO +A«bA)2B2 



A2 



KK' 



(11) 



Explicit expressions for the corresponding four-spinor 
eigenvectors | ^p± ) are too cumbersome to reproduce 
here, but will be used below to calculate the necessary 
matrix elements. In the case Akk' — ^ 0, | ip^ ) and | ) 



correspond to states adiabatically connected to | + |) 
and |— t) as B — >• 0, where -|-(— ) indicates r = 1 
(r = —1). More generally, the ± label in \tp±) simply 
indexes the two states which are degenerate as _B ^ 0. 

Substituting the expressions for \ip±) into Eq.© 
yields the rate (W^) for emission (absorption) of 
a phonon with momentum q and polarization a 



TTq 



hpLzOJq 
^2 



{N, + l)\{^+\Hl 



2^z dq 



27r ciLO 



nq 



hpL^UJ, 
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where Nq is the thermal occupation number for the 
phonon with hjjq = E_ — E^, and 



A; 



SO 



(Xq • s) T3 - florb (Xq • B) T3 



(7t 



(13) 

is the f7-transformed spin-phonon couphng operator of 
Eqs.® and ([6]) with the phonon operators and related 
dimensional constants removed. Using the exphcit ex- 
pressions for I ip± ) , we find 



(V+I/^.^IV--) 



2\Akk' 
A3A 



2i\A 



KK' 



A2 



M{q)AsOtJ'orbBz 



M{q) 



(14) 



where M{q) is the longitudinal form factor calculated in 
Appendix 

First, because A oc i?, it is apparent that both ma- 
trix elements in ([M]) vanish as i? — )- as expected for a 
transition between time-conjugate partners of a Kramers 
doublet. Second, both matrix elements vanish when 
Akk' — ^ 0, which is also expected because in this limit 
the two members of the Kramers doublet belong to op- 
posite valleys K and K', which are not coupled by long 
wavelength phonons. Finally, although it is not imme- 
diately obvious from Eq. p^ . in the limit B± the 
matrix elements for transitions involving the two polar- 
ization modes are equal in magnitude as required by 
axial symmetry: |( V+ \ \^P-)\ = \{^P+ \ \^-)\ ^ 

2\M{q)AsoAKK'f^orbB,\/A\ 

The relaxation rate F = 1 /Ti is the total rate of phonon 
emission and absorption in both channels x^: 



(15) 



To simplify the complicated expressions of Ea.(|14p. we 
take advantage of the fact that 



|Ms//Uort| < 1, 



(16) 



typically^Ti^ ~ 0.1/iorb, and that B± always appears 
multiplied by fiB, and so retain only the the leading term 
in fiorbBz- Using the energy conservation law huz = 
E- — E^ = hujq « 2fiorbBz and diu/dq ~ 2/3q, expression 
([T^ evaluates to 

j_^nmq)l 



2p 



AsqAkk' 
A2 



T 



3/2 



coth 



huz_ 
2fcsT 



(17) 

Not surprisingly, the uz dependence of expression pT)) 
has the same form as that found by Borysenko et al.— 
However, their relaxation rate F^g is suppressed by a 
very weak g-factor anisotropy Ag/g « 0.01—. Compar- 
ing prefactors, we find that the rates of relaxation due to 



these two mechanisms are related by the ratio 

2 



FAg 

F 



1 /Ag A^ 

2 ' AsoAkk' 



(18) 



For parameters in the regime which appears to be rele- 
vant for experiments [JHIi i-C- Aso ^ Akk', we find 
FAg/F ^ 1, indicating that the spin-orbit mediated 
mechanism studied herein should strongly dominate the 
behavior at low fields. 

We now provide an estimate of F based on Eg.lfTT]). 
with realistic experimental parameters. For /3, the coef- 
ficient of q^ in the flcxon dispersion, we use the result of 
Mariani and von Oppcni^ 



' 2fi{fi + A) 
{2fx + \)ps 



(19) 



where i? w i nm is ttie nanotube radm**, ps 
lO^'^g/cm^ is the graphene sheet density^!, related to 
the linear density p through p = 27ri?ps, and p ~ 
4A = 9 cVA" are Lame coefficientsi^. This gives 
hf3 sa 10** /leVnm^. Note that the theory of flexon vi- 
brations in nanotubes was first studied independently by 
Suzuura and Andoi^ and MahanSS^ who used different 
models of the curved graphene sheet. Mahan's theory 
and that of Suzuura and Ando, which was later applied 
by Bulaev et aii^ and by Mariani and von Opper^i^, 
predict different functional dependence of /3 on A and 
p. However, the numerical values obtained from each 
of these theories differ only by an insignificant factor of 
order 1. 

3/2 

In the low temperature limit, F is suppressed as uj^ 
'AS B ^ 0. The ratio of Ago to Akk' is strongly sample- 
dependent; for illustration we take Ago/ Akk' — 2. For 
a typical field B = lO'^ G nominally oriented along the 
tube axis, we find F w 300 s~^ for T ~ Q. For high tem- 
peratures the relaxation rate is enhanced by kBT/huz 
due to the presence of thermal phonons, and the ujz ^ 

1 /2 

limit of Eq. p7|) goes as F oc ■ 

In this section we focused on relaxation within the low- 
est Kramers doublet. When condition (fTBl) is satisfied, 
however, the results above can also be applied directly to 
relaxation within the upper Kramers doublet. 



IV. UPPER AVOIDED CROSSING 

We now turn our attention to spin relaxation in the 
vicinity of the upper avoided crossing circled in Figllji, 
where a minimum in Ti{B) was discovered by Churchill 
et ali^. Such an acceleration of spin-relaxation in the high 
temperature regime was first predicted in Ref. [l3 | for the 
case of electron-phonon coupling through the deforma- 
tion potential. Here we build on the idea outlined in 
that work, and consider how relaxation is affected in the 
presence of the deflection-coupling mechanism described 
above. As we will show below, for small level splittings 
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FIG. 2: (color online) Behavior of the spectrum in lim- 
iting regimes I^kk' = or B± — 0. a) Spectrum for 
AxK' ~ 0, B± 7^ 0. The levels Ei and E2 exhibit an avoided 
crossing with splitting controlled by the field misalignment 
B±. The crossing between E2 and E3 is exact, except when 
the field direction deviates very far from the tube axis (see 
text), b) Spectrum for Axj^i 7^ 0,B± = 0. Intervalley scat- 
tering opens an avoided crossing between levels Ei and E3. 
The exact crossing between Ei and E2 persists as long as 

^SO > 2{flB/tJ,orb)i^KK'- 



Lo the rate is parametrically enhanced by a factor of 
thus producing a robust enhancement of spin relaxation 
both in the high and low temperature regimes. 

Because the secular equation corresponding to Hamil- 
tonian ^ is fourth order, the spectrum of the system for 
moderate values of B cannot in general be described in 
any simple form. Furthermore, depending on whether 
or not the system possesses certain symmetries such 
as valley conservation {/S.kk' = 0) or axial symmetry 
{B±_ = 0), level crossings may be exact or avoided. The 
various possible regimes are reviewed in subsection IIV Al 
below (see also Ref.[l2j). In subsection IIV BI we calculate 
the relaxation rate close to the upper avoided crossing 
for the general case /^.k^i 7^ 0,i3^ ^ 0, and then pro- 
vide results for the special limiting cases ^kk' = and 
B^ = in subsection IIV CI 



A. Spectrum for B 



For B > 0, the spectrum of Hamiltonian (j?]) displays 
two avoided level crossings, shown in FigH^. Coupling 
between nominally degenerate states is provided by in- 
tervalley scattering A.kk' or by the magnetic field com- 
ponent perpendicular to the tube axis, B±. To provide 
a better understanding of what controls the splitting at 
each crossing individually, we first examine the two lim- 
iting cases Akk' = 0, B± ^ and Akk' 7^ 0, B± = 0, 
which arc described by simple analytical expressions. 

In the first case. Akk' ~ 0, B± 7^ 0, the four branches 



of the spectrum are given by 

El = fJ-orbB, 



E2 



i5q 



B, 



1 

2 
1 

2 
1 

2 

E4 — —HorbBz — - 

shown in Fig|5Ji. An exact crossing at the upper in- 
tersection can only occur if Ei = E2, which implies 
2^bBz = Aso,B± = 0; therefore the spectrum displays 
an avoided crossing for any B± 7^ 0. Similarly, an exact 
crossing at the lower (central) intersection occurs when 
E2 = E3, which requires 



(Aso 


- 2fiBB,)^ - 


\- AfilBl 


(Aso 


- 2/iBB,)2 - 




(Aso 


+ 2fiBB,)^- 




(Aso 


+ 2tiBB,y - 





cos^ ip „ — ^ sin^ (p 



^SO' 



(21) 



where (p is the angle between the magnetic field B 
and the tube axis z. This condition is always satis- 
fied for some value of B when the magnetic field di- 
rection does not deviate too far from the tube axis, 
tan^ tp < {norb/^J'Bf^ - 1 « 100. 

In the opposite case, Akk' 7^ 0, — 0, we find 

El = fiBB, +i 

E2 =-iiBB,+]^ 

E3 = ^^BBz — - 

Ea =—^bBz—- 



(Aso 


— 2florbBz 


2- 




(Aso 


+ 2florbBz 


2- 




(Aso 


— 2florbBz 


2- 




(Aso 


+ 2^lorbBz 


2- 


V^A\^, 



An exact crossing E\ = E2 exists at Bz satisfying 



4/z 



Tit 



A2 



A2 



(23) 



whenever 4^gA|.^, < - /i|)AgQ. Because /xb <C 

[lorbi this condition is unlikely to be violated in clean 
samples. The gap between E\ and E-^ persists for any 
Akk' ^ 0. 

More generally, both Akk' , B± 7^ 0, as in Fig[T^. Here 
we gain further insight by focusing on the regime 



2i-lbBz = Aso - Se, 



(24) 



where 6s is a small deviation from the upper cross- 
ing found in the Akk' , B± = limit. For small Se, 
Akk'i and Bx, we diagonalizc the Hamiltonian within 
the 2x2 subspace defined by the crossing levels, tak- 
ing into account coupling to the lower two levels to first 
order in perturbation theory. The perturbative treat- 
ment is justified by the large energy denominators com- 
pared to the coupling matrix elements, {^.orb^ ^■b)Bz S> 
6e, iibB±, Akk' ■ Because Akk' shifts the level crossing 
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point, which is otherwise determined by condition ([24)) . 
we also require ^kk' ^ (Morb/Ms )^so to ensure that 
the perturbation criteria are not violated. This latter 
restriction is rather mild, however, due to ^or& 3> /is- 

This treatment yields eigenvectors | ipi ) and | "02 ) (not 
shown) and eigenvalues 



£'1,2 — fJ-orbBz 



(25) 



where +(— ) is used for -Ei(2), sec Figl5^. Here 
+ 4(Mo.b-Ml)'B.VBSi. 



(26) 



In writing TZ, wc use Ea. (|24l) and keep only the lowest 
order term in 5e. In the limits IS.kk' = and/or B± = 0, 
Egs. ipSl) and (pS)) reduce to the earlier results ([201) and 
(El. 



B. Transition Rates for AjfK',B± 7^ 

Proceeding as in section Hill the amplitudes for transi- 
tions from I -01 ) to I ?/'2 ) due to coupling to phonons with 
polarization a are 



.Aso 



Miq), 



(27) 



A2 



with TZ as defined in Ea. (|26l) . In the expressions 
above, small corrections proportional to {B±/Bz)'^ and 
'/["iifJ-lrb ~ f^%)Bz] have been omitted. 
Interestingly, the contribution from y-polarized 
phonons is constant near the avoided crossing, with an 
amplitude equal to M(q)Aso/2. When B± = 0, the two 
amplitudes arc equal as required by axial symmetry; al- 
though the energy splitting Ei — E2 goes through zero in 
this limit, the amplitude remains finite. Unlike the case 
of the Kramers doublets here there is no time reversal 
symmetry to prevent such nonzero amplitude. 

Additionally, for ^ 0, the amplitude for emitting 
an x-polarized phonon vanishes for the detuning 5e = Jeq 
satisfying 



(28) 



Thus the coupling (and hence the relaxation rate) is 
asymmetric in detuning relative to the position of the 
anticrossing. For Akk' /^so ^ 1, the asymmetry is 
rather small. If it can be detected, however, it will allow 
the value of Akk' to be extracted. 

It is convenient to express TZ in terms of Seq, 

n = {fil,-fil)B,e{B,,B^), 



Here the energy splitting Ei~E2= HH^ylrb^ y-VjB^] 
e{Bz,B±) is controlled by a combination of B^ (de- 
tuning) and B± (avoided crossing gap). Summing the 
squares of the matrix elements (|27p . using the energy 
conservation law hujq = s{Bz, B±), and using the Golden 
Rule, Eq.dH), 



r(e) 



^so 



|M(g)|^ 



coth 



2fesr, 

(30) 

is the total rate of 



e(B„Bx) = J[5e{B,)-5eQY+AiilBl. (29) 



where T = Ea(W^a + W^) = ^/Ti 
transitions between | if^i ) and | "02 ) • 

As expected, the relaxation rate (|30p between states 
near the upper avoided crossing shows a resonant en- 
hancement proportional to e~^/^ (e^"^/^) in the low 
(high) temperature regime due to the combination of the 
divergent fiexon density of states at small energies and 
the deflection coupling mechanism. The singularity in 
the rate is a factor of 1/e stronger than that found by Bu- 
laev et al for the deformation potential coupling mech- 
anism. For larger values of e, the competition between 
the two mechanisms is sensitive to various system param- 
eters. Without an analytical expression for r(£) for the 
deformation coupling mechanism, however, it is difficult 
to make a more detailed comparison in this regime. 

Using the parameters of Ref.jH, Aso = 170 /ieV, 
e w 10 /icV, and T w 0.1 X, we find F « lO^s-i. De- 
spite the strong singularity for small energy splittings, 
this value is still 2 orders of magnitude smaller than the 
observed relaxation rate. However, the rate is highly sen- 
sitive to parameter values, including the fiexon tempera- 
ture and especially the avoided crossing splitting e. The 
latter depends on the angle of misalignment between the 
magnetic field and the tube axis, which was estimated to 
be 5° by electron micrograph, but if smaller could eas- 
ily lead to enhancement of the rate. Thus although the 
magnetic field dependence of Ti from this model supports 
the association of the Ti minimum observed by Churchill 
et ali^ with coupling to flexons, additional experiments 
are needed to understand the quantitative details of this 
relationship. 



C. Limiting cases Akk' = 0, Bx 7^ and 

Akk' / 0, Bx = 

To gain a fuller picture of the behavior near the up- 
per avoided crossing, we now examine the simple limits 
Aifif' = and B±^ = of section IIV Al In the limit 
Akk' = 0, B±^ 7^ 0, El and E2 exhibit an avoided cross- 
ing of splitting E1 — E2 ~ 2ij.bB±. At this point, Se = Seq 
and therefore ( V'2 | \ tpi ) vanishes; only flexons with 
y-polarization contribute to relaxation. The finite en- 
ergy splitting e at the avoided crossing, controlled by 
Bj_, provides a low-energy cutoff for the denominator in 
Eg. ([501) . Further away from the avoided crossing, the rate 
is suppressed by the decaying phonon density of states. 
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but gains an extra factor of 2 from the phonons with 
x-polarization that become active there. 

In the other Hmit, when I^kk' and = 0, there 
is an exact crossing when condition (f23| is satisfied. Here 
both flexon polarizations contribute equal transition ma- 
trix elements Aso/2. With constant matrix elements in 
the numerator and no low energy cutoff in the denomina- 
tor, expression (|30p diverges as B approaches the critical 
value defined by Eg. ([25]) . Of course, in any real experi- 
ment, the divergence will always be cut off by something, 
such as a small misalignment of the field axis. However, 
as found for the general case in the previous subsection, 
the behavior in all regimes displays an acceleration of 
spin-relaxation near the upper avoided crossing. This is 
a general property of the deflection coupling mechanism 
of spin-relaxation. 



V. DISCUSSION AND CONCLUSIONS 

Above, we found that the direct coupling between an 
electron's spin and the deflection of its host nanotube 
provides an efficient mechanism for electron spin relax- 
ation near narrow avoided level crossings where the en- 
ergy transfer is small. Efficient relaxation is made possi- 
ble by the diverging density of states for bending mode 
phonons with quadratic dispersioni^ii^ii^i^, oc g^. 
Such a dispersion law is predicted for clean systems such 
as perfect suspended nanotubes. What happens for dirty 
systems where the nanotube is placed on a substrate, 
and/or covered by an irregular coating left over from 
sample processing? 

The primary effects of disorder are twofold. Con- 
straints on the tube's motion alter the displacement pro- 
files of the normal modes. As a result, the spectrum of 
normal mode frequencies is altered, with spectral density 
at very low frequencies generally getting redistributed to 
higher frequencies. We begin this section by discussing 
normal mode profiles and providing some additional dis- 
cussion of the longitudinal form factor M{(\) that appears 
throughout the text. Then we discuss how changes in 
M{q) and dq/duj due to fiexon localization may affect 
the relaxation rates calculated above. 

By writing Eq.® for the nanotube deflection, we as- 
sume that the flexon normal modes are described by 
plane waves e*''^ . Even for a perfectly clean system, how- 
ever, the boundary conditions at the mechanically con- 
strained points of the tube can mix modes with momenta 
±g to yield normal mode proflles which are sums of sines 
and cosines. While these boundary conditions determine 
the nanotube length L^, the electronic length Ld is de- 
termined by a different set of constraints arising from 
electrostatic potentials created by external gates and/or 
impurities; thus it is not possible to find a completely gen- 
eral result for the form factor M{q), which depends on 
the shapes and relative displacements of both the elec- 
tron wavefunction and fiexon normal mode profile. In- 
deed the authors of Ref.[l9j concluded that asymmetric 




TT 27r Stt An 

qLd 



FIG. 3: (color online) Longitudinal form factor |M(g)p for 
symmetric square well confinement as described in Appendix 
\K\ with barrier separation Ld — 100 nm, for Eg = 0.03 eV 
(dashed blue) and Eg = 1 eV (solid red). For the large gap 
tube, kLd ~ TT and |M(g)p closely follows the hard wall limit, 
Eq. (|A16| ). with strong suppression for qLd > 2tt. For the 
small gap tube, softer confinement results in a larger effective 
length Ld, kLd = ^- When plotted against qLd, the profile 
nearly collapses onto the large gap result (dotted line). 



suppression of vibrational sidebands in transport through 
a nanotube quantum dot could be explained by the pres- 
ence of a vibrational mode localized near one end of the 
quantum dot. Thus we see that the general situation can 
be quite complicated; we consider the above-results for 
plane wave normal modes to be applicable for the "typi- 
cal," or "average" case. 

In section [TT] we pointed out that for small momentum 
(energy) transfer g — > 0, the factor M{q) I. For large 
momentum transfer qLd ^ 2tt, however, M{q) is strongly 
suppressed due to the fast oscillations of e*^^ relative to 
the electron wavefunction. Using the results of Appendix 
Elfor symmetric square well confinement, we plot |A'/((7)p 
for a large gap and a small gap tube in Figl31 For the 
small gap tube, where kLd ~ 0.75, confinement is softer 
and the electron wavefunction extends over a distance 
significantly larger than the barrier separation Ld- Thus 
suppression sets in at qL'^ 3> 2tt, where kL^ = tt. 

The asymptotic behavior of M{q) for large q contains 
two contributions. The classically allowed region pro- 
vides a contribution M{q) ^ 1/q'^ for qLd ^ 27r, see 
Ea. (jA16p . In the classically forbidden region, the expo- 
nential tails of the wave function contribute a Lorentzian 
decay law M{q) ^ l/{k'^ + q'^) with a prefactor that tends 
to zero as kcLd — >■ oo. Power law decay arises from the 
steps in the potential at z = ±Ld/2, but the suppression 
behavior is generic. For a smooth confinement potential 
one expects exponential decay, with a characteristic scale 
qL*d « 27r. 

The resonant enhancement of relaxation due to the 
flexon density of states singularity is thus not affected 
by M{q) as long as the minimum of the energy gap 
e{Bz,B±) = hfiq"^ at the upper avoided crossing corre- 
sponds to a phonon momentum q satisfying qL^ < 27r. 
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Using fi,/3 10^ /icV(nm)^ (sec Sec jllip and an energy 
gap of 10 ^eV as observed in Ref.[l|, we find qL'^ = 
1.6-3.2 for = 50-100 nm. Thus for the parame- 
ters of Ref. , the longitudinal overlap does not suppress 
relaxation at the upper avoided crossing, but will lead 
to suppression for larger energy splittings. Note that 
the same form factor appears in the matrix elements for 
the deformation potential couplingi^ and for the g-factor 
anisotropy deflection coupling^!, and thus does not affect 
the relative competition between these mechanisms. 

In the presence of disorder, the quadratically dispers- 
ing flexons are subject to the same conditions that lead 
to the familiar localization of all one-dimensional elec- 
tron eigenstates in a random potentia P^'^^ . For weak 
scattering characterized by ^ 1, where i is the lo- 
calization length, the density of states remains smooth 
and the normal mode profiles oscillate within envelopes 
with exponentially decaying tails, e~'^'/^. The exponen- 
tial envelope cuts off the rapid oscillations responsible for 
the suppression of M{q), allowing £ to substitute for Ld- 
Thus localization can in some cases help promote relax- 
ation at energy splittings larger than those where strong 
suppression by M{q) sets in for clean samples. Because 
the localization length £ is of the order of the mean free 
path22, which for short range scatterers is proportional 
to the phonon energy, £ (x q'^ and the clean limit (very 
weak scattering) should generally be recovered for large 
enough energy transfers. 

When scattering is strong, <C 1, the normal mode 
spectrum collapses to a collection of discrete resonances, 
and the normal mode profiles are highly distorted. In 
this case, the low energy modes responsible for the van 
Hove singularity that promotes spin relaxation near the 
upper avoided crossing can be destroyed. Rather than 
observing a smooth increase of the relaxation rate as the 
intersection is approached, a complicated non-monotonic 
system of resonances may then be found. 

It is also interesting to note that rigid rotational mo- 
tion of the entire sample, as described in FiglTjD, can be 
regarded as a global {q = 0) deflection of the nanotube. 
Consequently, simple laboratory vibrational noise can be 
a source of electron spin decoherence through the deflec- 
tion coupling mechanisni^. 

In summary, we have identified and analyzed an effi- 
cient mechanism of electron spin relaxation in nanotube 
quantum dots that results from spin-orbit-mediated cou- 
pling between the electron spin and nanotube deflec- 
tion. Due to the flexon density of states singularity at 
small energies, relaxation due to this deflection coupling 
mechanism is particularly efficient near level crossings 
or small avoided crossings. This mechanism is expected 
to dominate over other mechanisms such as deformation 
potentiali^ and anisotropic g-factor induced deflection 
coupling2i, which are suppressed by small phonon mo- 
mentum q and g-factor anisotropy, respectively. Finally, 
we predict a robust minimum of the spin relaxation time 
Ti near the upper avoided crossing in both the high and 
low temperature regimes, thus offering a flrm basis for 



understanding the observed Ti minimum of Ref. [sj . 

Note added in proof. Recently, electrical spin manip- 
ulation in curved carbon nanotubes based on a related 
spin orbit coupling mechanism was proposed^^. 
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Appendix A: Longitudinal Form Factor M{q) 

As mentioned in section |lll every matrix element 
( "0/ I i?s-ph I 'ipi ) includes a longitudinal form factor 

/oo 
dzn{z)e""', 
-oo 

which depends on the speciflc form of longitudinal con- 
finement for the quantum dot, but is independent of 
the composition of | ipiif) ) in terms of the basis states 
{| T .s)}. In this appendix we discuss the properties of 
M{q) in more detail, and provide an explicit result for 
M{q) in the case of square- well confinement. 

The property that a unique longitudinal overlap inte- 
gral M{q) can be factored out of all matrix elements of 
Hs-ph independent of the initial and final states relies 
on two key facts. First, the perturbation Hg-ph, Eq.(IS|), 
is diagonal in the isospin index. Thus for any initial 
and final states | tpi ) and \ipf)j which can be written as 
I V'n ) = X^T s Cn;rs| T s), whcrc n stands for i or /, 

( Ipf I Hs-ph I i^t ) =^C*f.^^Ci.r's' {t S\ Hs-ph \t' s') 
s.s 

= ^('*f:TsC^■,rs'{T s\H,-ph\T S').{A1) 

Furthermore, because each state |t s) = |t)(X)|s) is a 
product of orbital and spin parts, Eq.([5]) gives 

{ts\ if,_ph I r ) - ( r I e*«^ \t)-{s\ HJ_^^ (A2) 

where Hl_^^ is the spin-phonon coupling Hamiltonian of 
Eq.(IS]) projected onto the orbital state | r), with 

/oo 
dzoe*«^«(T|5(z-zo)|r) 
-oo 

factored out. Second, the density n(zo) = {t\5{z — 
zq ) I r ) , which includes a sum over densities on the two 
sublattices, is a real scalar. Because | r ) with r = 1 and 
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T = — 1 are time-reversal conjugate, n{zQ) is independent 
of T. Hence 



;^/|i?s-ph|^0 = A/(<7) J2 



ph 



(A3) 

thus completing the proof that a single, unique contribu- 
tion from the longitudinal degrees of freedom appears in 
all transition rates. 

In the presence of a magnetic field that breaks time 
reversal symmetry, the density nT,s(^o) ~ {t s\S{z — 
zo)\t s) is generally not independent of r and s. In 
particular, the term T^^orbB^ in Eq.(|3]) leads to a renor- 
malization of Eg of opposite sign in valleys K and K'. 
For the device of Ref . Q , with a gap Eg 30 meV and for 
fields up to 1.5 T, however, fiorbB/Eg < 1/30. Thus our 
calculation of the longitudinal wave functions to lowest 
order and the resulting factorization of M{q) is justified. 

We now calculate M{q) for a quantum dot in a car- 
bon nanotube formed by a symmetric step potential as 
shown in Fig 01 Although the details of M{q) depend on 
the form of the longitudinal confinement, this example il- 
lustrates its basic properties. Electron eigenfunctions in 
a square well potential were analyzed in detail in Ref. [l^l . 
Therefore here we focus on the main features and omit a 
step-by-step derivation. 

A pure "hard- wall" boundary condition cannot be de- 
fined for this system due to the relativistic dispersion 
and internal spinor structure (recall Klein phenomenon 
for relativistic particles in one dimension). The strongest 
confinement is realized when the potential step at the 
barrier is roughly equal to half the band gap, Vb ~ Eg/ 2. 
Here the tails of the wave function decay exponentially 
over a distance proportional to E~^ . For large gap tubes. 
Eg (X 1/i?, while for small gap tubes, the curvature- 
induced gap scales according to Eg oc cosirj/E?, where 
R is the tube radius and i] is the winding angle. 

The orbital states |t) with t = ±1 of Eq.([2]) corre- 
spond to the lowest quantized mode of longitudinal mo- 
tion in the K and K' valleys, respectively. Because we 
assume that the external potential V{z) only depends on 
the longitudinal coordinate z and not on the circumfer- 
ential coordinate x, these eigenfunctions of Hd can be 
factored as 



V't(z), 



(A4) 



where r = (cc, z), the circumferential wavevector is 
proportional to Eg, Eg = ifvupkc, and ipr is a two- 
component wavefunction describing amplitudes il!^{z) 
and V't (-2): oil the graphene A and B sublattices. Re- 
call that Kt = r27r/3a(l, for t = ±1 corresponds 
to the K {t ^ 1) OY K' {t ~ ~\) point of the graphene 
Brillouin zone, where a is the lattice constant. In terms 
of the pseudospin Pauli matrices cti and <T2 , 4'r satisfies 

hvp ( TCTifcc — icr2— ) + y{z) 



' dz 



^r = E^r, (A5) 



see Eq.([T]). In a homogeneous tube with constant V{z), 
the eigenstates are plane waves proportional to e**^^ 



V = Vo 



V = Vo 



Vo 


I 


II 


III 








V = 




E 


-id/2 


Ld/2 
Eg/2 = Vf-kc 



FIG. 4: (color online) We consider a quantum dot formed 
in a (small or large gap) semiconducting nanotube. The 
gate-induced potential V{z) is taken to be: V{z) = Vo for 
z < -Ld/2 OTZ> Ld/2, and V{z) = for -Ld/2 <z < Ld/2. 
The band gap Eg arises either due to quantization of trans- 
verse motion for pure semiconducting tubes, in which case 
Eg OC l/R, or due to curvature of the graphene sheet, in 
which case Eg oc 1/B? . Electron density n{z) is indicated by 
the red (dotted) curve. 



and the spectrum has a relativistic dispersion E = 

In the classically forbidden region I, z < —Ld/2, the 
potential takes the constant value V{z) = Vq. To empha- 
size the essential behavior and to simplify the equations, 
we pick the special value Vq = E, where E is determined 
self-consistently through the allowed wavevector k which 
satisfies the continuity relations at the boundaries, see 
Eq. ()A10[) below. With this choice, the wavefunction in 
the classically forbidden region is described by a single 
exponential living on only one sublattice: the spinor com- 
ponents V'^j and i/i^j satisfy 

<iW = ^(l-r)e'^^(^+^^/2), 
^U^) = ^(l + r)e^^(-+^''/2). (A6) 
Similarly, in region HI, x > Ld/2, 

Dr 



= ^(l + r)e-'==(^-^^/^), 
- ^(l-r)e-'^-(-^^/^). 



(A7) 



In these regions, the solutions in opposite valleys exist on 
opposite sublattices. However, the tota/ density "redzD is 
indeed the same for each case. 

In the classically allowed region II, V{z) — and the 
spinor components are plane waves ip^n oc e''^^re'^"^'= 
and V't II which satisfy: 



rkc — ik 
rkc + ik 



E 



(A8) 
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Here e^'P^ = - ik)/^/kfTl^, and E = Hvp y^kfTk^. 
Note that ip-k = —^k- 

Due to the degeneracy of states k and — fc, in region II 
a general state with energy E has the form: 



1 



— ikz ( 



1 



.(A9) 

Using the continuity condition of the wavefunction at 
the boundaries between different regions, we fix the rela- 
tive coefficients A^., BT,Cr, and Dt, and find the allowed 
wavevector k that gives a bound stated, 



tan kLri 



k 
he 



(AlO) 



which comes from the condition kLd = (pk mod tt. 

We calculate M{q) = Mi{q) + Afii(q) + Miii(g) piece- 
wise. In region I, Mi{q) — dz \Ar\^e 
evaluates to 



2p2fccZ,g(2fcc+i(j)z 



Mi{q) 



4Afsm^{kLd)e-' 
2kc + iq 



where = Ld/2 and the normalization constant sat- 
isfies 



1 4,2 
--: = ALd + -r- sm kLd - r 2fcid- 
J\ kc k 



In region III, 



Mm{q) 



AN sin'^ (kLd) e"'^ 



M*{q). 



(A12) 



(A13) 



2kc — iq 

Finally, in region II we must evaluate 

Mii{q)=AN J dx;[sin^fc(z + z,)+sin^fc(z-z,)] 6*"^"=, 

(A14) 



(All) which after some steps yields: 



Mu{q) = 4J\f 



2 sin qLd/ 2 



— cos kLd 



sm{k + q/2)Ld ^ sin(fc — q/2)Ld 



2k ^ 



2k- q 



r 



(A15) 



Because the density n{z) is symmetric in z, M{q) = 
Mi{q) + Mii{q) + Miu{q) is real, and even in q. 

For a high aspect ratio dot with kcLd ^ 1, we have 
kLd ~ TT. In this case, the contributions Mj and Mm 
from the tails vanish, and Eq. ljAlSp simplifies to 



Muiq) 



Stt^ sin(gLd/2) 
qLdi^TT^ - {qLdf 



(A16) 



which equals 1 for g = and decays as {qLd)~^ for qLd >• 
27r. This limit coincides with Eq.(56) of Ref . [l^ . Inter- 
estingly, for a small-gap nanotube with Eg w 30 meV 
and Ld = 100 nm as in Refs.jlQ, kLd ~ 0.757r, indicat- 
ing significant penetration into the classically forbidden 
region. 
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